Convergence in the QM‐only and QM/MM modeling of enzymatic reactions: A case study for acetylene hydratase

We report systematic quantum mechanics‐only (QM‐only) and QM/molecular mechanics (MM) calculations on an enzyme‐catalyzed reaction to assess the convergence behavior of QM‐only and QM/MM energies with respect to the size of the chosen QM region. The QM and MM parts are described by density functional theory (typically B3LYP/def2‐SVP) and the CHARMM force field, respectively. Extending our previous work on acetylene hydratase with QM regions up to 157 atoms (Liao and Thiel, J. Chem. Theory Comput. 2012, 8, 3793), we performed QM/MM geometry optimizations with a QM region M4 composed of 408 atoms, as well as further QM/MM single‐point calculations with even larger QM regions up to 657 atoms. A charge deletion analysis was conducted for the previously used QM/MM model ( M3a , with a QM region of 157 atoms) to identify all MM residues with strong electrostatic contributions to the reaction energetics (typically more than 2 kcal/mol), which were then included in M4 . QM/MM calculations with this large QM region M4 lead to the same overall mechanism as the previous QM/MM calculations with M3a , but there are some variations in the relative energies of the stationary points, with a mean absolute deviation (MAD) of 2.7 kcal/mol. The energies of the two relevant transition states are close to each other at all levels applied (typically within 2 kcal/mol), with the first (second) one being rate‐limiting in the QM/MM calculations with M3a ( M4 ). QM‐only gas‐phase calculations give a very similar energy profile for QM region M4 (MAD of 1.7 kcal/mol), contrary to the situation for M3a where we had previously found significant discrepancies between the QM‐only and QM/MM results (MAD of 7.9 kcal/mol). Extension of the QM region beyond M4 up to M7 (657 atoms) leads to only rather small variations in the relative energies from single‐point QM‐only and QM/MM calculations (MAD typically about 1–2 kcal/mol). In the case of acetylene hydratase, a model with 408 QM atoms thus seems sufficient to achieve convergence in the computed relative energies to within 1–2 kcal/mol.Copyright © 2013 Wiley Periodicals, Inc.     

Biologically Active 1-Aminodeoxy and 1-O-Alkyl Derivatives of The Powerful D-Glucosidase Inhibitor 2,5-Dideoxy-2,5-Imino-D-Mannitol

ABSTRACT

By an Amadori rearrangement of easily available 5-azido-5-deoxy-D-glucofuranose with dibenzylamine and subsequent catalytic hydrogenation of the resulting 5-azido-1-(N,N-dibenzyl)amino-1,5-dideoxy-D-fructopyranose, 1-amino-1,2,5-trideoxy-2,5-imino-D-mannitol was obtained in only two steps and in excellent overall yield. Likewise, other amines were employed to introduce extended side chains ultimately suitable for attachment of the inhibitor to solid supports. The reported rearrangement reaction is a high yielding, convenient and apparently general entry to 1-aminodeoxyketopyranoses modified at C-5, facilitated by the ring enlargement of the aldofuranose to the ketopyranose as an additional driving force. A range of selected chain extended analogues was prepared by acylation of N-1. Inhibitors obtained exhibit K _i-values with D-glucosidases in the micromolar range. Interestingly, 1-N-acylation resulted in superior inhibitory activities, as did the addition of a hexyl chain.     

Nonlinear wave propagation in locally dissipative metamaterials via Hamiltonian perturbation approach

Nonlinear Dynamics - The cellular microstructure of periodic architected materials can be enriched by local intracellular mechanisms providing innovative distributed functionalities. Specifically,...     

Oligosaccharides: Defense Inducers,Their Recognition in Plants,Commercial Uses and Perspectives

Plants have innate immune systems or defense mechanisms that respond to the attack of pathogenic microorganisms. Unlike mammals, they lack mobile defense cells, so defense processes depend on autonomous cellular events with a broad repertoire of recognition to detect pathogens, which compensates for the lack of an adaptive immune system. These defense mechanisms remain inactive or latent until they are activated after exposure or contact with inducing agents, or after the application of the inductor; they remain inactive only until they are affected by a pathogen or challenged by an elicitor from the same. Resistance induction represents a focus of interest, as it promotes the activation of plant defense mechanisms, reducing the use of chemical synthesis pesticides, an alternative that has even led to the generation of new commercial products with high efficiency, stability and lower environmental impact, which increase productivity by reducing not only losses but also increasing plant growth. Considering the above, the objective of this review is to address the issue of resistance induction with a focus on the potential of the use of oligosaccharides in agriculture, how they are recognized by plants, how they can be used for commercial products and perspectives.     

Taba: A Tool to Analyze the Binding Affinity

Evaluation of ligand-binding affinity using the atomic coordinates of a protein-ligand complex is a challenge from the computational point of view. The availability of crystallographic structures of complexes with binding affinity data opens the possibility to create machine-learning models targeted to a specific protein system. Here, we describe a new methodology that combines a mass-spring system approach with supervised machine-learning techniques to predict the binding affinity of protein-ligand complexes. The combination of these techniques allows exploring the scoring function space, generating a model targeted to a protein system of interest. The new model shows superior predictive performance when compared with classical scoring functions implemented in the programs Molegro Virtual Docker, AutoDock4, and AutoDock Vina. We implemented this methodology in a new program named Taba. Taba is implemented in Python and available to download under the GNU license at https://github.com/azevedolab/taba . © 2019 Wiley Periodicals, Inc.     

Competition between Liquid-liquid De-mixing,Crystallization, and Glass Transition in Solutions of PLA of Different Stereochemistry and DEET

Liquid-liquid(L-L)de-mixing and vitrification of solutions of either crystallizable poly(L-lactic acid)(PLLA)or non-crystallizable poly(D/L-lactic acid)(PDLLA)with 50 m%N,N-diethyl-3-methylbenzamide(DEET)were analyzed by calorimetry and cloud-point measurements,which allows drawing conclusions about the effect of polymer stereochemistry on the phase behavior.Regardless of the PLA stereochemistry,vitrification of the solutions on fast cooling,hindering crystallization of PLLA,occurred below-20℃ and suppressed prior L-L de-mixing.The experimental results prove that crystallization in samples containing crystallizable PLLA,observed at around 55℃ on slow cooling,is not preceded by L-L de-mixing.     

pH Tuning of Water-Soluble Arylazopyrazole Photoswitches

Arylazopyrazoles are an emerging class of photoswitches with redshifted switching wavelength, high photostationary states, long thermal half-lives and facile synthetic access. Understanding pathways for a simple modulation of the thermal half-lives, while keeping other parameters of interest constant, is an important aspect for out-of-equilibrium systems design and applications. Here, it is demonstrated that the thermal half-life of a water-soluble PEG-tethered arylazo-bis(o-methylated)pyrazole (AAP) can be tuned by more than five orders of magnitude using simple pH adjustment, which is beyond the tunability of azobenzenes. The mechanism of thermal relaxation is investigated by thorough spectroscopic analyses and density functional theory (DFT) calculations. Finally, the concepts of a tunable half-life are transferred from the molecular scale to the material scale. Based on the photochromic characteristics of E- and Z-AAP, transient information storage is showcased in form of light-written patterns inside films cast from different pH, which in turn leads to different times of storage. With respect to prospective precisely tunable materials and time-programmed out-of-equilibrium systems, an externally tunable half-life is likely advantageous over changing the entire system by the replacement of the photoswitch.     

Aqueous Diels–Alder reactions for thermochemical storage and heat transfer fluids identified using density functional theory

Thermal storage and transfer fluids have important applications in industrial, transportation, and domestic settings. Current thermal fluids have relatively low specific heats, often significantly below that of water. However, by introducing a thermochemical reaction to a base fluid, it is possible to enhance the fluid's thermal properties. In this work, density functional theory (DFT) is used to screen Diels–Alder reactions for use in aqueous thermal fluids. From an initial set of 52 reactions, four are identified with moderate aqueous solubility and predicted turning temperature near the liquid region of water. These reactions are selectively modified through 60 total functional group substitutions to produce novel reactions with improved solubility and thermal properties. Among the reactions generated by functional group substitution, seven have promising predicted thermal properties, significantly improving specific heat (by as much as 30.5%) and energy storage density (by as much as 4.9%) compared to pure water.